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ABSTRACT 


In this work, a gradual solar energetic particle (SEP) event observed by multi¬ 
spacecraft has been simulated. The time profiles of SEP fluxes accelerated by an 
interplanetary shock in the three-dimensional interplanetary space are obtained by 
solving numerically the Fokker-Planck focused transport equation. The interplan¬ 
etary shock is modeled as a moving source of energetic particles. By fitting the 
1979/03/01 SEP fluxes observed by Helios 1, Helios 2, and IMP 8 with our sim¬ 
ulations, we obtain the best parameters for the shock acceleration efficiency model. 
And we also find that the particle perpendicular diffusion coefficient with the level of 
~ 1% - 3% of parallel diffusion coefficient at 1 AU should be included. The reservoir 
phenomenon is reproduced in the simulations, and the longitudinal gradient of 
SEP fluxes in the decay phase, which is observed by three spacecraft at different 
locations, is more sensitive to the shock acceleration efficiency parameters than that is 
to the perpendicular diffusion coefficient. 

Subject headings: Sun: activity — Sun: coronal mass ejections (CMEs) — Sun: particle 
emission 
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1. INTRODUCTION 


There are two categories of solar energetic particle (SEP) events: impulsive events and 
gradual events. The impulsive events are small, last for hours, are rich in electrons, 3 He and heavy 
ions, have relatively high charge states, and are produced by solar flares. In contrast, gradual 
events are large, last for days, are electron poor, have relatively low charge states, and are related 
to the shocks driven by interplanetary coronal mass ejections (ICMEs). Some large ICME driven 
shocks can cover a large range of solar longitudes and latitudes in the interplanetary space, and 
the observers located at different locations can be connected to different parts of the shocks by 
interplanetary magnetic field (IMF). Therefore, multi-spacecraft observations by, e.g., Helios 1 


and 2, usually show a large difference of SEP time profiles at different longitudes ( Re a me s et ah 


1996, 


19971) . In addition, in some SEP events, both of solar flare and shock can exist in the 


same event. Therefore, there are also shock-flare-mixed SEP events which can be identified 
by abundance ratio and charge state of heavier ions. 

Multi-spacecraft observation data provide essential information to understand the processes 
of particle acceleration and transport in the heliosphere. Multi-spacecraft observations in the 
ecliptic plane, e.g., by Helios 1 and 2, or at different latitudes and radial distances, e.g., by 
ACE and Ulysses , usually show two interesting phenomena in some gradual events. For the first 
one, SEP events could be observed by multi-spacecraft with a very wide spatial distribution that 


could be wider than the size of the source (e.g., 


Wiedenbeck et al. 


20131) . This phenomenon can 


be explained in part as a result of the effec t of particle perpendicular d i ffusion, and has bee n 


investigated in detail with simula tions (Zhan g et al. 


Wang et al. 


2012; 

Oin et al. 

2013; 

Droge et al. 

2014; 


2009!; 


He et al. 


2011 


Dresing et al. 


Gomez-Herrero et al 


2015 


2012; 


Wang & Oin 


20151) . For the second one, the SEP fluxes observed by widely separated multi-spacecraft 


usua 

19971 


l y show similar intensit i es within a small ~ 2-3 factor in different positions (Reames etah 


McKibben et al. 


2001 


Lario et al. 


2003 


Tan et al. 


2009; 


Zhang et ah 


2009; 


Oin et al. 
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19921) named it as “reservoir”. Recently, 


2013; W ang & Oi nl |20 1 51) . iMcKibbenl (119721) discovered thi s phenomenon, and iRoelof et al. 


Wang et al. 


(2012) use a numerical code (denoted as 


Shock Partic 
the model of 


e Transport Code, S PTC) con s iderin g shock as a moving SEP source by adopting 


Kalle nrode ( 19971) . lOin et al. 


(20131) reproduced the reservoir phenomenon wit 


different shock acceleration efficiency and perpendicular diffusion. In addition, 


Wang & Oin 


( 20151) investigated the spatial and temporal invariance in the spectra of gradual SEP events. In 


their simulations, SPTC is used with the IMF set as Parker field model, and the disturbance of 
the IMF caused by ICME is ignored. They found that shock acceleration efficiency, parallel 
diffusion, adiabatic cooling, and perpendicular diffusion are four important factors in forming the 
reservoir phenomenon, and the first three factors are the main factors with the last factor being a 
secondary one. The peaks of SEP fluxes are mainly controlled by shock acceleration efficiency 
and parallel diffusion. And the fluxes decay in the similar ratio due to the effect of adiabatic 
cooling. Furthermore, because of the effect of perpendicular diffusion, the longitudinal gradient in 
the SEP fluxes, which is observed by spacecraft located at different locations, is further reduced. 


Observationally, the four factors change significantly in different SEP events (Kallenrode 


1996 


1997), so that only in the gradual SEP events when the values of the controlling effect parameters 


are appropriate can the reservoir phenomenon be formed. 

Generally, SEP acceleration by shocks are calculated in two major approaches: in the first 


a pp roach, SEPs are injected at t 


1995 


Ka: 


Oin et al 


(Lee 


lenrode & Wibberenz 


1983: 


re shock with an assumed injection 


1997 


Lario et al. 


1998 


Kallenrode 


2001 . 


stren gth (He ras et al. 


Wang et al. 


1992, 


2012; 


2013 ), while in the second approach the acceleration o f SEPs by shocks are included 


Gordon et al. 


1999; 


ZuoetaL 


2011 


ZuoetaL 


20131) . Each approach has its own 


advantages, the first one can provide a reasonable description of the SEP fluxes in gradual 
events by focusing on the transport of energetic particles without the thorough knowledge of the 
details of shock acceleration, while the second one can help us to better understand diffusive 
shock acceleration by including more physics details. There are also some studies attempting 
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to combine the advantages of the two approaches. The SEPs acce 


processes 


Rice et al, 


are included in these studies (Ng et al. 


1999 


Zank et al. 


eratio n and transpo rt 


200 3). 


2000; 


Li et al. 


2003; 


Because the calculation of SEP flux needs a precise mechanism for particle’s injection into the 
diffusive shock acceleration, which is currently not available, we adopt the first approach to inject 


SEPs at the shock with an assumed i njection strength, so that we could 


shock accelerated particles’ transport (Wang et al. 


2012; 


Oin et al. 


2013 


ocus on interplanetary 


Wang & Oin 


20151) . 


In this paper, as a continuation of our previous research, we further study the the shock 
acceleration efficiency and SEPs transport in a gradual SEP event by comparing our 
simulations with multi-spacecraft observations. We describe the SEP transport model and 
the shock model in section 2. We show the observation results in section 3. We show the shock 
geometry and the effect of perpendicular diffusion in Section 4. We show the simulation results 
and their comparison with multi-spacecraft observations in Section 5. We summary our results in 
Section 6. 


2. MODEL 


In this work we follow previous research < 

C 

>in et al. 2 

SD 
O 

o 

Oil; Oin et a 

1. 201 

3; 

Zhang et al. 

2009 

Droge et al 

2010 

; Zuo et al. 

2011; 

Zuo et al 

2013; 

Wang et al. 

2012 


Wang et al. 

2014) 

to model the transport of SEPs. One can write a three-dimensional focused 

transport equation as (Skilling 

1971; 

Schlickeisei 

2002 

; Oin et al 

. 2006 

Zhang et al. 

2009 

) 


df 

dt 
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where f(x,p,p,t ) is the distribution function averaged over gyrophase, t is the time, x is the 
position in a non-rotating heliographic coordinate system, p, v, and p are the particle momentum, 

A 

speed, and pitch-angle cosine, respectively, in the solar wind frame, b is a unit vector along the 
local magnetic field, V sw = V sw r is the solar wind velocity, and the magnetic focusing length L 
is obtained from the divergence of the IMF background strength B 0 , i.e., L - [b -VlnBoj • Here, 
we use the Parker field model for the IMF. The equation CD) includes almost all important particle 
transport effects, such as particle streaming along the field line, adiabatic cooling in the expanding 
solar wind, magnetic focusing in the diverging IMF, and the parallel and perpendicular diffusion 
coefficients. 


We define a pitch a ngle diffusion coefficient following (IBeeck & Wibberenzl Il986l: 


Oin et al. 

2005; 

Oin et al. 

2006) 


D mm (p) = D 0 v(R L k min y- 2 (p s ~' + h) (l - p 2 ) , 
where the constant D 0 is adopted from Teufe l & Schlickeiser (2 00311 


( 2 ) 


r* (SB s iab\ n(s- i), ^ 

Do ~( Bo ) As ^ min ’ (3) 

here 6B s i ab /B 0 is the magnetic turbulence level of slab component, K, = pc / (\q\ B {) ) is the 

maximum particle Larmor radius, q is the particle charge, l s i ab is the slab turbulence correlation 

length, k mm = 1 /l s [ ab is the lower limit of wave number of the inertial range in the slab turbulence 

power spectrum, and s = 5/3 is the Kolmogorov spectral index of the magnetic field turbulence in 

the inertial range. The constant h arises from the non-linear effect of magnetic turbulence on the 


pitch angle diffusion at p = 0 (Oin & Shalchi 
h = 0.01, and k m ; n = 32 AU -1 . 


2009 . 


2014). In the following simulations, we set 


Following I Jokipiil (ll966l) . lHasselmannl (11968). and 


Earl (1974), the relationship between D 


MM 


and parallel mean free path (MFP) dy is written as 


3v r +l 

~8 J-i 


(1 ~P 2 f 


D, 


dp. 


MM 


( 4 ) 































-7- 


In addition, the parallel diffusion coefficient Ky is related to Ty by /ey = vTy/3. 

The perpend icular diffusion coefficient is set by using the nonlinear guiding center theory 


(Matthaeus et al. 

2003 

) with the following analytical approximation ( 

Shalchi et al. 

2004, 

2010) 


1 

3 V 


So ) ^ 2s r(f + i) 


,2/3 


h d 


A, 


1/3 


I~bb 


(5) 


where B 2D /B Q are the turbulence level of 2D component, and l 2D is the correlation length. T 
is the gamma function. I is a unit tensor. In our simulations, l 2D is set to 3.1 x 10~ 3 AU, 
(6B 2D ) 2 / (SB s i ab ) 2 = 4, and 5 = 5/3. As a result, the values of parallel and perpendicular diffusion 
coefficients can be altered by changing the magnetic turbulence level 6B/B 0 = ^SB 2 lab + 6B 2 D /B 0 . 

The particle injection into the sho ck at position (r, 0, tp) and time t with mom entum 


p is specified by 


boundary values from 


Wang et al. 

(2012 

), and 

Oin et al. 

(2013) 


Kallenrode & Wi bbe renz (119971) : 


Kallenrode (200lb : 


f b (r, e,(p,p,t) 
S(r,8,ip) 

m<p) 


a ■ 6(r - v s t ) • S ( r , 6, ip) ■ p y • ^(0, ip) 


r 

10(0, ip) 1 

1 exp 

0c 


[ 1 if \<f>(6, p)\ < <p s 

0 otherwise, 


( 6 ) 


where r is the solar radial distance, v s is the shock speed, v s t = r b0 + n ■ Ar with n = 0,1,2, • • •, n 0 
and A r being space interval between two ‘fresh’ injections, r b0 is the inner boundary, S is the 
shock acceleration efficiency which specifies the particles ejection and changes with a power law 
in radial distance and is exponential towards the flank of shock, r c is the radial normalization 
parameter, a and cf) c are the shock acceleration efficiency parameters, determines the spatial scale 
of shock front, (f> is the angle between the center of shock and the point at the shock front where 
the particles injected, and <p s is the half width of the shock, y is the energy spectral index of SEP 


source. 
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In order to numerically solve the transport eguation ([]]). a time-bac 


cward Markov stochastic 


process method is used by following Zh ang ( 1999b , see also 


Oin et al. 


(20061) for details of the 


application of the methods to study SEPs. Here, we denote the numerical code to calculate the 
propagation of energetic particles, which treats the CME driven shock as a moving particle source, 
as Shock Particle Transport Code, i.e., SPTC. 


3. OBSERVATIONS 

The SEP data from the Helios University of Kiel particle instrument provides proton 
and heavier nuclei with energy range between 1.3 and several hundred MeV/nuc, and 
electrons of 0.3 - 0.8 MeV. The SEP data from the IMP 8 Goddard Medium Energy (GME) 
Experiment provides an energy range of 0.5 - 450 MeV proton and 2 - 450 MeV/nuc heavier 
nuclei, and relativistic electrons. During the time period March 1-11, 1979, a gradual SEP 
event was observed by Helios 1 and 2, and IMP 8, which were located near 1 AU ecliptic, but 
at different longitudes. On March 1, 1979, the radial distances of Helios 1, Helios 2, and 
IMP 8 on March 1,1979 are 0.95 AU, 0.93 AU, and 0.99 AU, respectively. And the Helios 1, 
Helios 2, and IMP 8 were located at 17°, 57°, and 84° in the heliographic inertial coordinate, 
respectively. Figures [DEI and[3]show the time series of SEP fluxes, interplanetary magnetic field, 
and solar wind measurements from the Helios 1, Helios 2, and IMP 8, respectively. In each 
of the three figures, from top to bottom, the observation data are for SEP fluxes, magnetic field 
strength |B|, polar and azimuthal field angles 6 and (p in the Selenocentric Solar Ecliptic (SSE) 
coordinate for Helios 1 and 2 but in the Geocentric Solar Ecliptic (GSE) coordinate for IMP 8, 
plasma density N, plasma temperature T, and bulk solar wind speed V. The vertical line indicates 
a interplanetary shock passage. From Figured] we find that during the period March 3-5 Helios 
1 observed both a shock and an ICME event characterized by increases in the solar wind speed, 
density, temperature, and the rotation of the azimuthal field angle <p. From Figured] we find that 
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Helios 2 only observed a shock without ICME. The ICME driven shock was detected by the 
Helios 1 and Helios 2 during March 3. If the same shock was detected by IMP 8, it should 
be detected nearly in the same time interval. However, during March 3, we find that IMP 8 
observed neither a shock nor an ICME in Figure |3j This SEP event is also investigated by 


Lario et al. (2006) and 


Reames (2010). In these two studies, they also found that IMP 8 did 


not observe a shock. 


Just before the onset of the March 1-11, 1979 gradual SEP event at 1 AU, a solar flare is 


observed (Kallenrode et al, 


1992 ). The flare is located at S23E58, and the time of the maximum of 


soft X-ray is at 10:19 on March 1. If we assume the time of the maximum of soft X-ray to be the 
moment of formation of the CME driven shock, the average speed of the shock can be obtained, 
i.e., 0.57 AU/day and 0.47 AU/day for the average speed of shocks detected by Helios 1 and 2, 
respectively. In this work, we further take an average of 0.57 AU/day and 0.47 AU/day, and set the 
shock speed to be a constant 0.52 AU/day. 


4. SHOCK’S GEOMETRY AND PERPENDICULAR DIFFUSION 

In fact, the geometry of an interplanetary shock would be very complicate. For 
simplicity, we do not include shock angle (obliquity) in our model, but only use a cone to 
model the shock’s geometry in our simulations. The shock nose (the center of the cone) is 
located in ecliptic, and the longitude of the shock nose is set as the same as that of solar flare. 
In this event, the shock nose points to E58. Since the locations of spacecraft are known, 
then we can determine the relative positions between the shock nose and the spacecraft. In 
this case, the Helios 1, Helios 2, and IMP 8 are 9° east, 31° west, and 58° west of the shock 
nose, respectively. Although the direction of the shock nose is a hypothesis, we can use the 
spacecraft’s in-situ observation to identify if the hypothesis is reasonable. In this SEP event, 
Helios 1 detected the shock and ICME, but IMP 8 didn’t. These observations implied that 
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the Helios 1 is located the nearest to the shock nose, and IMP 8 is the farthest. Therefore, 
the locations of spacecraft in our model are consistent with the in-situ observations of 
spacecraft. Furthermore, the solid angle of shock front is set as 35°, which agrees with 
the fact that both Helios 1 and Helios 2 detected shock, but IMP 8 did not. According to 
the observations of Helios 1, there is no significant difference in the SEP fluxes when the 
spacecraft is inside and outside the ICME. As a result, the disturbances of IMF caused by 
ICME could be ignored, and the IMF is set as Parker spiral in our simulations. 

Based on our shock model and the observations of Helios 1, Helios 2 , and IMP 8, we plot 
a cartoon for illustrating the cross-section of the shock and the locations of three spacecraft in 
Figure [4j In this figure, the shock front is indicated by the dashed arc line, and the shock nose is 
indicated by the dashed-arrow radial line passing through the center of the shock. The big solid 
circles indicate the locations of the three observers, Helios 1 and 2, and IMP 8, and the small 
ones indicate the protons of SEPs. As the shock propagates outward, the shock front is connected 
with the Helios 1 by the IMF first, then with the Helios 2, at last with the IMP 8. Note that the 
spacecraft’s field line is not connected to the shock front all the time. For example, the field line of 
the IMP 8 is not connected to the shock at the beginning, then as the shock front moves to a larger 
radial distance than that of the IMP 8, the field line becomes connected to the shock. Eventually, 
the observer will be disconnected from the shock as the shock continues to propagate outward. 

In Figure [4] (a), we assume that SEPs propagate in the interplanetary space without 
perpendicular diffusion. In this case, the SEPs can only propagate along the IMF. The particles can 
be detected at the onset time when the spacecraft’s IMF is connected to the shock front. In Figure 
|4](b), however, we assume that SEPs propagate in the interplanetary space with perpendicular 
diffusion. In this case, SEPs can cross magnetic field lines. The particles can be detected before 
the spacecraft is connected to the shock by field lines, and they can more easily spread in the 
interplanetary space. In March 1, 1979 SEP event, the onset time of SEP fluxes observed by three 
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spacecraft was very close. Therefore, in order to reproduce the observations with our model, 
perpendicular diffusion must be included in the following simulations. 


5. THE RESULTS OF SEP SIMULATIONS AND THEIR COMPARISONS WITH 

OBSERVATIONS 

The parameters used in the simulations are listed in tables Q] and [2] The parameters in table 
(Hare the same in all of following simulations, but the ones in tableware different in individual 
simulations. In the observations, the speed of solar wind is always changing with time, but 
the variance in the speed would not change our main conclusions significantly. As result, we 
used a constant speed of 400 km/s in our simulations. 


5.1. Effect of Perpendicular Diffusion on SEP Flux 

Figure [5] shows the simulation and observation results of the time profiles of SEP fluxes. 
The black and red lines indicate the observations of 3 - 6 MeV protons and the simulations of 
5 MeV protons, respectively. And the solid, dotted, and dashed lines are corresponding to the 
observations of Helios 1, Helios 2, and IMP 8, respectively. The vertical lines indicate the 
moment when the spacecraft’s field lines are connected to the shock. Due to the different locations 
of the three spacecraft and the width of the shock, the spacecraft’s field lines are connected to 
the shock front at different time. The Helios 1 is connected to the shock by IMF first, then is the 
Helios 2, at last is the IMP 8. The simulations in Figure [5] (a) are noted as Case 1. Note that the 
parameters for different cases are shown in Table El The onset time when the shock is connected 
to IMP 8 by IMF is nearly two days later than that to Helios 1. With perpendicular diffusion, 
particles can cross magnetic field lines in the interplanetary space, so they can be detected before 
the spacecraft’s magnetic field line is connected to the shock. As a result, the onset time of SEP 
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fluxes in the simulations is very close because of the effect of perpendicular diffusion. As we 
can see, the simulations roughly agree with the observations in Figure [5] (a). The simulations in 
Figure [5] (b) are noted as Case 2. Figure [5] (b) is the same as Figure [5] (a) except that the ratio of 
perpendicular diffusion coefficient to parallel one is larger. In Figure 0(b), SEP flux observed by 
Helios 1 is very close to that in Figures [5] (a). However, during the early stage of flux rising phase, 
SEP fluxes observed by Helios 2 and IMP 8 in Figure [5] (b) are slightly larger than that in Figure 
|5](a). This is because all particles arrive at Helios 2 and IMP 8 by crossing field lines before the 
observers field lines are connected to the shock front. Due to the effect of stronger perpendicular 
diffusion in Figure [5] (b) than that in Figure [5j (a), during the decay phases, the spatial gradients 
of SEP fluxes in Figure [5] (b) are slightly smaller than that in Figure [5] (a). In this SEP event, the 
effect of perpendicular diffusion is very important when the spacecraft’s field line is disconnected 
from the shock front, and it also helps to reduce the spatial gradient in the decay phase of SEP 
fluxes. 


5.2. Effect of Source Injection Profile on SEP Flux 

Figures [6] (a) and (b) are the same as Figure [5] (a) but with different shock acceleration 
efficiency parameters, and the simulations in Figures [3 (a) and (b) are noted as Case 3 and Case 
4, respectively. In simulations of Figure [6] (a), the shock acceleration parameter a is set to 2.5. 
Comparing with the simulations in Figure [5] (a), the shock acceleration efficiency decreases more 
slowly with radial distance, and the longitudinal gradient in the simulation fluxes is larger. In the 
simulations of Figure [6] (b), the shock acceleration parameter <p c is set to 5°. Comparing with the 
simulations in Figure [5] (a), the shock acceleration efficiency decreases more quickly toward the 
shock flank, and the simulation fluxes increase more quickly during the rising phase. The peak 
intensity of SEP fluxes are mainly determined by shock acceleration efficiency, and the gradient 
of SEP fluxes in the decay phase is more sensitive to the shock acceleration efficiency parameters 
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than the perpendicular diffusion coefficient. In this SEP event, at the peak time of flux for Helios 
2 (IMP 8), the flux for Helios 2 (IMP 8) was similar to that for Helios 1 which had been in decay 
phase. The reservoir phenomenon is formed in the decay phase in the SEP fluxes. 


5.3. Effect of Adiabatic Cooling on SEP Flux 

Figures [7] is the same as Figure [5] (a) but with different injection spectrum. The simulations 
in Figure 0 is noted as Case 5. When propagating in the heliosphere, SEPs would lose energy 
because of the effect of adiabatic cooling, so that the SEPs with the same energy observed by 
spacecraft later generally originate with higher energies at the source. Because the source has 
a negative energy spectral index y, the flux decreases more quickly with a smaller y. The fluxes 
observed by Helios 1, Helios 2, and IMP 8 decay as a similar ratio due to the effect of adiabatic 
cooling. In the simulations of Figure [7] the injection spectral index y is set to -8.5. Comparing 
with the simulations in Figure [5] (a), the simulation fluxes decrease more quickly, and the decay 
ratios of SEP fluxes are mainly determined by the effect of adiabatic cooling and the energy 
spectral index y. 


5.4. Higher Energetic Particles 

In Figures [H the black and red lines indicate the observations of 32 - 45 MeV protons 
and the simulations of 35 MeV protons, respectively. The simulations in Figure [8] is noted as 
Case 6. In simulations of Figure [8l the shock acceleration parameter a is set to 4. Comparing 
with the simulations in Figure |5] (a), the peak time of simulation fluxes in Figures [8] comes 
earlier. This means that source acceleration efficiency decreases with radial distance faster 
at higher energy channel. 
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6. DISCUSSION AND CONCLUSIONS 

In this work, we study a gradual SEP event which was observed by Helios 1 and 2, and IMP 
8. The event studied in this paper began on March 1, 1979, and lasted nearly eight days. By 
solving a three-dimensional focused transport equation, the fluxes observed by the three spacecraft 
are calculated. The transport equation we use in this work includes many important particle 
transport effects, such as particle streaming along the field line, solar wind convection, adiabatic 
cooling, magnetic focusing, and the diffusion coefficients parallel and perpendicular to the IMF. 
By comparing the simulations and the observations, we get the shock efficient parameters and 
diffusion coefficients to get the best fitting. The following are our major findings. 

The effect of perpendicular diffusion is very important when the spacecraft’s field line is 
disconnected with the shock front. In 1979/03/01 SEP event, the in-situ observation shows that an 
ICME is detected by Helios 1, but not by Helios 2 or IMP 8 which are in the west of Helios 1. 
Therefore, Helios 1 is located near the center of shock front but Helios 2 and IMP 8 are located 
in the west flank of the shock, if it is assumed that the ICME is located behind the center of the 
shock. Furthermore, the interplanetary shock is only observed by Helios 1 and 2, but not by IMP 
8. In the case without perpendicular diffusion, the IMP 8 can not detect any SEP when the radial 
distance of the shock is smaller than 1 AU. However, according to the observations of IMP 8, the 
particles can be detected by IMP 8 before IMP 8’s field line is connected to the shock front. In 
order to reproduce the observations of 3 - 6 MeV and 32 - 45 MeV protons in the simulations, the 
perpendicular diffusion coefficient with the level of about 1% - 3% of parallel one at 1 AU should 
be included, 

The peaks of fluxes are mainly determined by the shock acceleration efficiency, and the 
decay ratios of SEP fluxes are mainly due to the effect of adiabatic cooling and the energy spectral 
index. Due to the effect of adiabatic cooling, the fluxes observed by Helios 1, Helios 2, and IMP 
8 decay as a similar ratio. The shock acceleration efficiency S is assumed as r~ a exp (- \(p\ /<p c ) in 
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the simulations. The shock acceleration efficiency is weaker in the flank of the shock front than 
that in the center, and decreases as the solar radial distance r increases. For 3-6 MeV protons, 
by comparing the simulations and the 1979/03/01 SEP event observed by Helios 1, Helios 2, 
and IMP 8, we find that a is set to 3, <p c is set to 10°, and the energy spectral index y is set to 
-6.5. We also calculated the SEP fluxes in the cases of different model parameters which are 
listed in the Table |2j For 32 - 45 MeV protons, we find that a is set to about 4 in the Case 
6 of Table |2j Comparing with the simulations in the Case 1 of 3 - 6 MeV protons, the a is 
larger in the case of 32 - 45 MeV protons, which means the source acceleration efficiency 
decreases with radial distance faster at higher energy channel. 

The longitudinal gradient of SEP fluxes in the decay phase observed by different spacecraft 
is more sensitive to the shock acceleration efficiency parameters than the perpendicular diffusion 
coefficient. In 1979/03/01 event, at the peak time of flux for Helios 2 ( IMP 8), the flux for 
Helios 2 {IMP 8) was similar to that for Helios 1 which was in decay phase, and the reservoir 
phenomenon was formed in the decay phase of SEP fluxes. In this event, the reservoir 
phenomenon is also reproduced in the simulations. However, in some other SEP events, at the 
peak time of flux for Helios 2 {IMP 8), the flux for Helios 2 (IMP 8) was significantly different 
than that for Helios 1 which was in decay phase. The reservoir phenomenon can not be formed in 
the decay phase of SEP fluxes with normal diffusion coefficients ( k ± /k\\ is a few percent). 

The IMF is set to the Parker field model in our simulations, and the disturbance of the IMF 
caused by ICME is ignored. However, according to the observations of Helios 1 in this event, 
there is no obvious difference in the SEP fluxes after the spacecraft enter the ICME comparing to 
the fluxes detected by Helios 1 when it is out of ICME. As a result, we assume that the disturbance 
of the IMF would not change the main results of our simulations. Furthermore, our model 
can not be used to study the shock-flare-mixing SEP events. In the future work, we plan to 
include the particle source combined with shock and flare to model the shock-flare-mixed 
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SEP events. And we also intend to include a realistic three-dimensional ICME shock with the 
disturbance of IMF caused by ICME, so that the SEP acceleration and transport in the heliosphere 
can be investigated more precisely. 
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Helios 1 



1979 Mar 


Fig. 1.— Solar wind plasma and magnetic field parameters observed at the Helios 1 spacecraft. 
From top to bottom, magnetic field strength |B|, polar and azimuthal field angles 6 and 0 (in SSE 
coordinate), plasma density N, plasma temperature T, and bulk solar wind speed V. The vertical 
line indicates a interplanetary shock passage. 
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Helios 2 



Fig. 2.— Same as Figure [T] except that the parameters are observed at the Helios 2 spacecraft. 
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IMP-8 






3-6 MeV 
6-11 MeV 
11-23 MeV 

29-43 MeV 

43-51 MeV 


Fig. 3.— Same as Figure Q] except that the parameters are observed at the IMP 8 spacecraft and 
olar and azimuthal field angles are in GSE coordinate. 























































































































-24- 


Without perpendicular diffusion With perpendicular diffusion 



Fig. 4.— Geometry of a shock cross-section with three spacecraft at different locations. The 
dashed-arrow radial line indicates the center of the shock. The big solid circles represent the three 
spacecraft, and the small ones represent protons of SEPs. 
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Fig. 5.— Comparison of the observations of 3 - 6 MeV proton fluxes with black lines and the 
simulation results of 5 MeV protons with red lines. The observers are located at 1 AU in the 
ecliptic, but at different longitudes. The black lines show the time profiles of the observation 
fluxes, and the red lines indicate the simulation fluxes. The vertical lines indicate the moment of 
the spacecraft’s field line is connected to the shock. 
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64 66 68 

Doy in 1979 


Fig. 6.— Same as Figure [5] (a) but with different shock acceleration efficiency parameters a and 

0c- 



Fig. 7.— Same as Figure [ 5 ] (a) but with different injection spectrum y. 
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Fig. 8.— Comparison of the observations of 32 - 45 MeV proton fluxes with black lines and 
the simulation results of 35 MeV protons with red lines. The observers are located at 1 AU in 
the ecliptic, but at different longitudes. The black lines show the time profiles of the observation 
fluxes, and the red lines indicate the simulation fluxes. 


Table 1: Model Parameters Used in the Calculations. 


Parameter 

Physical meaning 

Value 

ysw 

solar wind speed 

400 km/s 

r c 

radial normalization parameter 

0.05 AU 

A r 

shock space interval between two fresh injections 

0.001 AU 

Vs 

shock speed 

0.52 AU/day 

(f>s 

shock width 

35° 

VbO 

inner boundary 

0.05 AU 

n,\ 

outer boundary 

50 AU 
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Table 2: Model Parameters Used in the Calculations. 


Case 

a 

0c 

V 

K ± a 

r 

0 SB/B 0 ) 2 

1 

3.0 

10° 

0.48 AU 

1.2% X AT|| 

-6.5 

0.4 

2 

3.0 

10° 

0.24 AU 

3.1% X AT|| 

-6.5 

0.8 

3 

2.5 

10° 

0.48 AU 

1.2% X AT|| 

-6.5 

0.4 

4 

3.0 

5° 

0.48 AU 

1.2% X AT|| 

-6.5 

0.4 

5 

3.0 

10° 

0.48 AU 

1.2% X AT|| 

-8.5 

0.4 

6 

4.0 

10° 

0.66 AU 

0.97% x /C|| 

-6.5 

0.4 


for protons in the ecliptic at 1 AU. 



